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Abstract 

A new approach for for estimating the acoustic pressure in the near field of a mi- 
crophone based on non-intrusive direct measurement of acoustic particle velocity is 
proposed. 

This method enables the estimation of the acoustic pressure inside a domain 
located in front of the microphone membrane. The acoustic pressure is calculated 
using the acoustic particle velocity on the frontiers of this domain and a physical 
model based on the Green function of the system. 

Results are obtained using the acoustic velocity measured with Particle Image 
Velocimetry (PIV) in front of a microphone excited with a plane wave inside a 
rectangular waveguide. They show that the diffraction of the plane wave by the 
microphone leads to an increase of the acoustic pressure on the microphone edge in 
the order of magnitude of 0.1 dB. 

Key words: PIV, pressure estimation, acoustic particle velocity measurement, 
integral formulation. 
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1 Introduction 



Over the years, several methods have been developed for microphone calibra- 
tion (for a review see [1]). These methods can be classified into two sorts : (i) 
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relative calibration which provides an estimate of the sensitivity of a micro- 
phone as a function of a reference sensitivity and (ii) absolute calibration [2] 
which leads to an estimation of the sensitivity without any reference micro- 
phone. 

For absolute calibration, the reciprocity technique is usually used which pro- 
vides a typical precision of around 0.05 dB in an enclosed field configuration. 
In a free field configuration, this technique has been adapted [3,4] and stan- 
dardized [5] . Nevertheless, free field absolute calibration suffers from numerous 
problems which are not resolved at present: the location of the acoustic center 
of the microphone is crucially important [6] , the generation of standing waves 
between the different microphones perturbs the measurements and, generally, 
external reflections disturb the acoustic field in the vicinity of the microphone. 

Microphone calibration using non-intrusive optical techniques such as Laser 
Doppler Velocimetry is proposed by some authors. These studies are mainly 
conducted in enclosed field conditions, more precisely in waveguides excited 
with plane waves. Two approaches are used. 

On the one hand, authors consider that the impedance of the medium is 
known and measure the acoustic velocity at a single point. This supposes that 
the boundary conditions of the medium are perfectly known. We call it the 
"global approach". In this case, the acoustic pressure can be estimated at 
the measurement point thanks to the impedance. Acoustic pressure can also 
be estimated elsewhere using a propagation model of the system under study. 
Taylor [7], MacGillivray [8,9] and Koukoulas et al [10] use this approach. First 
results are encouraging: Taylor [7] shows that this method is accurate within 
±0.03 dB at 500 Hz, Mac Gillivray et al [8,9] reach an accuracy of around 
0.1 dB and Koukoulas et al [10] propose an accuracy of 0.2 dB at 170 Hz. 

On the other hand, authors consider a volume of fluid and measure the acoustic 
velocity on the volume boundaries. A propagating model of the fluid enables 
the estimation of the acoustic pressure everywhere in the volume. We call this 
the "local approach". For plane waves a slice of fluid is considered and the 
acoustic velocity measured at two points to estimate the acoustic pressure in 
the middle of the slice. Degroot et al [11] use this approach and show that using 
a (u — u) probe with two LDV measurements provides a minimum uncertainty 
on the pressure estimation of 0.013 dB for frequencies of 1360 Hz and 680 Hz. 

The approaches described above for enclosed field calibration could be used 
in a free field. The first (global) approach only requires the measurement of 
the acoustic velocity at a single point but also requires that the impedance 
of the fluid for free field conditions is known. This can only be the case if 
the boundary conditions of the system are perfectly known, for example in a 
semi-infinite domain which can be reproduced with a semi anechoic chamber. 
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The second (local) approach does not require the boundary conditions of the 
system under study to be known but instead the acoustic velocity must be 
measured at a number of locations near the microphone membrane, in order 
to estimate the acoustic pressure on the membrane. 

The local approach requires the characterization of the acoustic nearfield of a 
structure by LDV or other non-intrusive techniques. Previous work has been 
done by Gazengel et al [12], who measure the acoustic particle velocity in 
front of a loudspeaker to characterize its acoustic radiation in a free field 
using the LDV technique. Schedin et al [13] propose measuring the acoustic 
particle velocity in the vicinity of a plate using two-reference-beam double- 
pulsed holographic interferometry and Moreau et al characterize the acoustic 
field in a waveguide boundary layer using LDV and PIV techniques [14]. 

In this paper, we propose a preliminary study for calibrating microphones in 
free field conditions using a non-intrusive measurement technique and the lo- 
cal approach described above. Section two presents the general method used 
for estimating the acoustic pressure on the microphone membrane using the 
acoustic velocity measured at many locations near the membrane. Section 
three presents the experimental study of the acoustic field in the close vicinity 
of a microphone. In this section, the acoustic particle velocity is measured 
using a PIV technique. In the fourth section, a 2D model of the acoustic field 
near the microphone in derived using the Green function of the volume under 
study. Finally, the acoustic pressure field is estimated in front of the micro- 
phone membrane using the measured acoustic velocity and some discussion is 
given. 



2 General formulation of the acoustic pressure field in the vicinity 
of a microphone membrane 

In this section, a microphone excited with a plane wave is studied. An analyt- 
ical development of the pressure field is proposed, considering a fluid domain 
(air) located in front of the microphone membrane (see Fig. 1). Acoustic pres- 
sure inside this volume can be estimated at the position (r, 9, z) = (w, z) by 
means of measurement of normal velocity on the volume boundaries with an 
integral formulation. 

The system under study (see Fig. 1) is an air fluid column with length L and 
a circular section of radius r a . The circular sections (Si, S2) and the surface of 
the fluid column S 3 are subjected respectively to the normal acoustic velocities 
v ni , v n2 and v n3 . As described in Fig. 1 and according to experimental results 
presented in section 3, the normal acoustic velocities v ni and chosen to 

be uniform on the circular section and v n:j is considered to vary on the section 
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Fig. 1. Studied volume in front of the microphone. 



A general integral formulation of the acoustic pressure in the volume V is 
proposed hereafter. 

In the frequency domain the acoustic pressure p(f, t) is written as 

p(r,t)=p(r)eP ut , (1) 

where u> is the acoustic pulsation defined by u> = 2nf. In the linear acoustic 
approximation and in the case of a perfect fluid, the acoustic propagation in 
the air fluid column V, is given by 

(A + k2)p(r) = 0, (2) 

where k = oj/cq is the wave number and Co is the sound celerity. The boundary 
conditions associated with eq. (2) are written 



= -jupv ni (r) for r 6 (0, r a ), 9 e (0, 2vr), z = 0, (3) 
-jupv n2 {r) for r G (0,r o ), 9 G (0,2vr), z = L, (4) 
-jupv nz (r) for r = r a , 9 G (0, 2n), z G (0, L), (5) 



dn 
dp(r) 

dn 



where p is the air density and d/dn = d n is the normal derivative to the 
surface S. The acoustic pressure field p(f) in the volume V (described by the 
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closed surface S) is written, in this integral form, as [] 

p(r) = JJJ G(f,f )f(f )dV + J J [G(r,r )dn p(r ) - p(f )dn G(f,f )}dS , 

(6) 

where the function f(fo) describes the sources distributed inside the volume 
V and G(f,f ) is the Green function defined by the follonwing equation 

(A + fc2)G(r, f ) = -5(f- f ) in V, (7) 

and the boundary conditions on the surface S. In eq. (7), 8(f — r*o) is the Dirac 
distribution. The boundary conditions for the Green function are chosen as 



d no G(r,r ) = for r £ (0,r a ), 9 £ (0,2vr), z = 0, (8) 
d no G{r, fb) = for r £ (0, r a ), £ (0, 2vr), z = L, (9) 
9 no G(f,r ) = for r = r a , 9 £ (0,2tt), z £ (0,L), (10) 

to allow a description of the acoustic pressure field as a function of the acoustic 
velocity on the boundaries. Using eqs. (7), (8), (9) and (10), eq. (6), without 
any sources (f(r ) = 0), is written 

p(r) = IJ s G(r,r )dn oP (r )dS . (11) 

Eq. (11) combined with eqs. (3), (4) and (5) leads to the following expression 
for the acoustic pressure in the volume V 

where i = 1, 2 and 3. Knowing the acoustic velocity on the surface Si, S2 and 
S3 and the Green function of the system, the acoustic pressure can determined 
in the volume V with the eq. (12). 



3 Experimental characterization of the acoustic velocity field in 
the vicinity of a microphone membrane 



In this section, the acoustic velocity field is estimated experimentally using 
Particle Image Velocimetry (PIV) inside the domain close to the microphone 
membrane as shown in figure 1. The microphone is positioned inside a waveg- 
uide. PIV has previously been used to study acoustic flow in a waveguide [14]. 
The acoustic velocity field is measured inside a plane region (laser sheet) of 
dimensions 2 r a , L. 
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3. 1 Experimental set-up 



3.1.1 P IV system 



The PIV system uses a pulsed Cu-Laser Oxford Lasers LS20 — 50 with a time 
between pulses of 20 /is. The mean power of the laser is 20 W and each pulse 
has a duration of 25 — 30 ns with a wavelength of 510.6 nm (green). The 
beam from the laser is converted to a light sheet for delivery to the region of 
interest. A CCD camera (Sensicam Double Shutter) is used for the acquisition 
of the PIV images, with a resolution 1280 x 1024 pixels. A Berkley Nucleonics 
model 500 A pulse generator is used in conjunction with National Instruments 
Lab VIEW software to allow the capture of PIV image pairs at any phase in the 
period of the acoustic cycle. 30 measurements are made at a given phase and 
averaged to estimate the velocity field. Cross-correlation and post-processing 
are carried out on in-house PIV software. 



3.1.2 Acoustic system 



The experimental set-up is made up of a JBL 244QH loudspeaker mounted on 
a closed Perspex tube (with 10 mm wall thickness) with a length L = 0.5 m 
and a square section S — 0.1 x 0.1 m 2 . The first cut-off frequency is 1720 Hz 
(first transverse mode). In the waveguide, a 1 inch B&K microphone is placed 
parallel to the guide axis. Fig. (2) shows the acoustic set-up and the PIV 
system. 



A stationary plane wave with a frequency / = 680 Hz is established in the 
waveguide. The microphone membrane is located between a node and an antin- 
ode of acoustic velocity to measure a sufficient acoustic velocity amplitude 
for PIV. The light sheet position is adjusted so as to graze the microphone 
membrane, in a plane corresponding to the diameter of the membrane. The 
waveguide is "seeded" using a SAFEX fog machine. The seeding is introduced 
during 3 to 5 seconds and 10 to 15 minutes are required before doing the mea- 
surement. The CCD Camera is positioned perpendicular to the light sheet and 
focused on the illuminated fog particles. The observation window corresponds 
to a rectangular section of the volume under study (dimensions 2r a x L). Mea- 
surements are realized for 20 equally spaced phase steps in the acoustic period 
(see Fig. 3). 
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Fig. 2. Experimental set-up 



0.05 
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Fig. 3. View of the 20 measurement phases during one acoustic period. 
3.2 Experimental results 

3.2.1 General analysis 

The acoustic velocity field measured for the phase 3 of the acoustic period 
is shown on the Fig. 4. The z-axis is the guide axis. The dimensions of the 
PIV image are 0.02 x 0.024 m 2 . This PIV image highlights the evolution of the 
acoustic velocity field streamlines at the vicinity of the microphone membrane. 
In this representation, the velocity vectors go toward the membrane for a 
positive velocity amplitude. The shape of the field lines are approximately 
symmetrical around the microphone axis. 
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Fig. 4. View of experimental acoustic velocity field in the vicinity of the microphone 
for phase 3. 

3.2.2 Longitudinal acoustic velocity 

Fig. 5 shows the longitudinal acoustic velocity function of r for z = m 

and for z = 0.015 m. In this figure, the microphone membrane is located at 
z = 0.015 m. This figure shows that the longitudinal velocity amplitude is 
almost constant for 2 = 0. The acoustic wave can therefore be considered 
to be plane at this position (corresponding to surface Si, see Fig. 1) in the 
domain under study. However, velocity v n \ shows a variation of ±5.10 -3 m.s -1 
around the mean value (5.10 -2 m.s -1 ), which introduces an uncertainty in the 
pressure estimation. 

In the very closed vicinity of the microphone (the region defined by r £ 
[—0.004; 0.004] m and z G [0.013; 0.015] m), the acoustic velocity field vanishes 
due to the membrane stiffness. For r = 0, the amplitude of the acoustic veloc- 
ity v z decreases from 0.045 m.s -1 {z — m) to 2.10 -3 m.s -1 (z = 0.015 m). 
Assuming the membrane velocity equals the acoustic velocity at z = 0.015 m, 
this result shows that the membrane velocity is very small compared with the 
acoustic velocity (ratio ~ 1/22) measured at z = (15 mm from the mem- 
brane). This result should be confirmed by complementary direct measure- 
ments of the membrane velocity. In further works, the microphone membrane 
velocity should be measured by means of a Laser Vibrometer in order to con- 
firm this hypothesis. 
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Fig. 5. Acoustic velocity amplitude 1) _2 clS £L function of r for z = m (a) and 
z = 0.015 m (b). 

3.2.3 Radial acoustic velocity 

Figure 6A shows the radial acoustic velocity amplitude 1J r clS db function of r 
for z = m and z = 0.015 m. In the very close vicinity of the membrane, 
the acoustic velocity amplitude v r increases with r. At z — 0.015 m and for 
r = m, the radial acoustic velocity amplitude is v r = 1.10~ 4 m.s -1 and for 
r = 0.011 m, v r = 0.45 m.s -1 . v r is maximum for r = 0.012 m which illustrates 
the acoustic leakage at the edge of the microphone. In the vicinity of the 
microphone, the distribution of the radial velocity amplitude is not exactly 
symetric (Fig. 6A). This asymetry can be due to microphone misalignment, 
error in the velocity estimation with PIV or a weak air current caused by 
thermal effects. For greater distance from the microphone (z < 0.015m), the 
radial velocity amplitude becomes symetric as shown on Fig. 6B. 

Finally, Fig. 7 illustrates the radial acoustic velocity amplitude ClS db func- 
tion of z for r = m and r = 0.012 m corresponding to surface S3 shown in 
Fig. 1. The radial acoustic velocity amplitude increases from zero (r = m) 
to 0.033 m.s -1 (r = 0.012 m). This phenomenon illustrates again the presence 
of acoustic leakage on the edge of the microphone and has to be taken into 
account in equation 12. 

Taking these experimental results into account, we consider in the following 
that the longitudinal velocity v ni is constant, that the membrane velocity v n2 
can be neglected (v n2 = 0) and that the radial velocity v n3 depends on the z 
coordinate. The profile of function of z is estimated analytically in 

the following (§5). 
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Fig. 6. (A) Radial acoustic velocity amplitude function of r for z = m 

(a) and z = 0.015 (b). (B) Radial acoustic velocity amplitude function of r 

for (a) z = 0.0136 m (...), (b) z = 0.0139 m ( ), (c) z = 0.0142 m (*) and (d) 

z = 0.0145 (+) 




4 Analytical model of the acoustic pressure field in the vicinity of 
a microphone membrane 



In this section, we derive a specific expression of acoustic pressure in volume 
V using the normal velocity profiles v nx and v m estimated from experimental 
results (§3). We assume an axial symmetry of the system, which allows the 
derivation of a 2D model of the acoustic field. 
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4-1 Green function of the system 



The Green function G(r,r ) solution of the problem described by eq. (7), can 
be presented as 

(Atf + 5 Z + k2)G(w, z\ w , z ) = —5(w - w )5(z - z ) in V, (13) 

where = (l/r)d r (rd r ) + (l/r2)dg2. The boundary conditions (eqs. (8), (9) 
and (10)) are expressed as following 



d no G(r, f ) = d Zo G(w, z\ w , z ) = for z = 0, (14) 
d no G(r, f ) = d Z0 G(w, z; w , z ) = for z = L, (15) 
d no G(r, f ) = d ro G(w, z; w , z ) = for r = r a .. (16) 

The Green function G(w, z; Wq, Zq) can be written as a discrete sum of eigen- 
functions ^^(w) of the cylinder V under the form 

oo 

G(w, z; w , zq) = g»v(z,zo)^»v{wo)^»v{w), (17) 

H,u=0 

where g^ v {z, Zq) are dependent on the position z and ^^(w) are solutions of 
the following problem 

(A* + k wlxv 2)-$ ^{w) = -5{w - wo) Vr G (0, r a ) and V0 G (0, 2vr), (18) 
drj&^w) = for r = r a and V# G (0, 2tt). (19) 

The eigenfunctions ^^(w) = ^^(r, #) take the following form 

*^(r, 0) = V cos(//0) J^{k wllv r), (20) 

where J M is the Bessel function and the eigenvalues k w ^ v are given by 

k Wilv = ^, with J'^ u ) = 0. (21) 

'"a 

The coefficients A^ u are found using the orthogonality of the eigenfunctions 
and are expressed as 

A (22) 

^ (1 + ^)^2(1-^)72^(7^)' 

Using eqs. (13), (17) and (18), the coefficients g^ u {z,zo) are solutions of the 
following relation 

(5 Z 2 + k2 xiu ,)g iw (z, zq) = 6(z - zq) for z G (0, L) (23) 
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with k2 g/ja , = k2 — k Wfll/ 2. The boundary conditions are given by 

d z g^{z, zo) = for z = 0, (24) 
9 z g^{z, zo) = for z = L, (25) 

and the solution is written 



cos(k zllu z)cos[k zllu (z - L)] 
9^{z, z ) = 1 . , ^ n if z < zo, (26) 

cos(k ZfMl ,zo) coa[k Zfll/ (z - L)] 
9nv{z, zo) = —r — 7- if z > Zq. (27) 



4-2 General formulation of acoustic pressure in the volume V 



The acoustic pressure in the volume V can now be determined by means of the 
Green function of the system and the normal acoustic velocity on the surface 
S. Using eq. (17) in eq. (12), the pressure is written as 

~ ~ oo 

P(r) = -jupY jj Y g»u{z,zo)^^{wo)^»v{w)v ni dSi, (28) 

i Si /i,u=0 

for % = 1, 2 and 3 corresponding respectively to the surface Si and the normal 
velocities v nr Setting 



Pi(r,9,z) = -jujp Y // 9^{z,Zii)^^(wo)^^{w)v ni dSi (29) 



the total pressure in the volume V is written as 

3 

p(r,6,z)=YMr,9,z) (30) 

i=l 

where % indicates the considered surface of the volume. The pressure due to 
the different surfaces can now be calculated separately to show the influence 
of each surface area. 



4-2.1 Calculation of the pressure field pi(r, 8, z) 

Using eqs. (20) and (21) in eq. (29), the pressure field pi(r,9,z) is expressed 
as 
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Pi(r,9,z) =-jup g IJ ,v(z,0)A2 tlv cos(ii9)J ll (—r)2'K6, 



mO 



H,v=0 



v ni (ro)J^(— r Q )r dr d9 (31) 
o r a 

where v ni (ro), the normal acoustic velocity field for x = 0, is independant of 
9. Taking into account the cylindrical symmetry of the system, the acoustic 
pressure field p\ is also independent of 9 (which implies that p = 0) and is 
written as 



pi(r, z) = -jup g 0u (z,0)A2 0u J (— r)2ir [ v ni (r )Jo(—ro)rodr d9o, 

r a Jo r a 

(32) 



u=0 

with 



cos(k z0u (z — L)) 

9ou{z, 0) = — — — for z > 0, (33) 

k z0u sm(k z0u L) 

and 

A2 0v = I ■ (34) 

r a 2TvJ ('jQ U ) 



4-2.2 Calculation of the pressure field p2{r, 9, z) 

The acoustic pressure field p 2 (r,9,z), considered as independent of 9 (due to 
the cylindrical symmetry), can be expressed as 



P2(r,z) = -jupY] g 0u (z, L)A2 0u J (—r)27r f v n2 (r )J (—r )r dr d9 , 

r a Jo r a 

(35) 



u=0 

with 



cos(k z0u z) 

g 0u (z, L) = - for z < L, (36) 

K zOv Sm[K z0l/ L/) 



and A2 0u defined by eq. (34). 



4-2.3 Calculation of the pressure field p%(r, 9, z) 

The normal acoustic velocity v n3 on the surface S3, considered as independent 
of r and 9 (due to the cylindrical symmetry), depends only on the z coordinate. 
The acoustic pressure field p 3 (r,9,z), considered as independent of 9 (due to 
the cylindrical symmetry), is expressed as a function of the normal acoustic 
velocity v n3 (z) on the surface S3 in using eqs. (20) and (21) 

00 lOu f L 

Pa(r,z) = -jujp V A2q u J ('Jou)Jo(— r)2nr a / v n3 (z )gou(z, z )dz . (37) 

t^O r a J 
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The expression for go„(z, z ) depends on the acoustic velocity profile v n3 (z) on 
the surface S3. 



4-3 Calculation of acoustic pressure in the volume V 



Using eqs. (30), (32), (35) and (37) and, as shown by the experimental results 
(section 3.2), considering a incident plane wave on the surface Si and v n2 (r) = 
on the surface 5*2, the total pressure field in the volume V can be written as 



. . . cos[fc(z - L)] 2up ( ~ Jo(^r) 1 



k sin kL 1 r a J (7o^) k z0l/ sm(k z0u L) 

cos[k zQu (z - L)\ / v n3 (z Q ) cos[k z0u z Q }dz 
Jo 



L 



+ cos[k Z Q U z\ J v n3 (z ) cos[k z0u (z - L)]dz 1 J , (38) 

where A 0u 2 and go„(z, z ) have been substituted by their expressions given by 
eqs. (34), (24) and (25). 

To express the total pressure field in the volume V, the acoustic velocity profile 
v n3 (z) on the surface 5*3 must be determined. The experimental results (section 
3.2) suggest that this profile can be modeled by a parabolic curve defined as 

v nz (z) = -f]z2 for z G [0, L] and 77 > 0. (39) 



Finally, the calculation of the integral function in eq. (38) using eq. (39) leads 
to the following result for the total pressure in the volume V (see appendix 
??) 



. . . cos[H,2 — L)} 
P ^ Z)= JUP k sin kL V 

x <^ z2sm(k z0l/ L) 



.2up 



00 j o( ^ r) 



E 

\u=Q 



n 



Jo(n/ov) k2 z0v sm(k z0v L) 



sin(k z0 „L) + Lcos(k z0u z) 



<*zOi> 



40) 



where k = k 



zOO- 



Now, the characterization of the acoustic pressure field in the vicinity of the 
microphone (volume V) can be made with the knowledge of the normal acous- 
tic velocities on the volume boundaries. 
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5 Results and discussion 



The calculation the acoustic pressure in front of the microphone requires the 
determination of the coefficient of the parabolic curve describing the variation 
of the normal acoustic velocity on the surface S3. This coefficient 77 defined in 
eq. (39) is determined by a minimization method. Fig. 8 shows a comparison 
between the experimental data and the model (defined by eq. (39)) with 77 = 
135 m~ 1 .s~ 1 . PIV measurements of the acoustic velocity amplitude normal to 
5*3 have been obtained for the phase 3 of the acoustic period (see Fig. 3). 

0.035 1 . 1 1 




Fig. 8. Normal acoustic velocity amplitude on the surface S3 as a function of z, 
experiment (— ) and model (• • • ) defined by eq. (39). 

Figs. 9 and 10 show the calculation of the pressure amplitude in the volume 
V for v = 2 modes and for v = 90 modes respectively [y > corresponds 
to radial modes defined in eq. 40). The influence of the number of modes on 
the calculation of the acoustic pressure is clearly shown by comparing these 
two figures. The diffraction effects on the acoustic wave due to the micro- 
phone are visible even with 2 modes but this phenomenon is described more 
precisely when the modes number increases. According to the experimental 
results, the presence of the microphone distorts the velocity field streamlines 
allowing an acoustic leak on the border of the microphone (Fig. 4). This phe- 
nomenon leads to an increase of the acoustic pressure amplitude on the edge 
of the microphone. The acoustic pressure amplitude is not constant along the 
microphone membrane, varying from 126.65 dB at the microphone center to 
126.8 dB at the microphone border. Due to the volume dimension and to the 
acoustic frequency (/ = 680 Hz), only the plane mode is propagative and 
all the higher modes are evanescent. The result of the calculation shows that 
the microphone edge creates diffraction of the plane wave which leads to the 
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redistribution of the acoustic energy of all the modes and consequently to the 
excitation of higher modes. 




Fig. 9. Calculation of acoustic pressure amplitude in the vicinity of the microphone 
using eq. (40) with v = 2. 




-0.012 -0.003 -0.004 0.00^ 0.008 0.012 

Fig. 10. Calculation of acoustic pressure amplitude in the vicinity of the microphone 
using eq. (40) with v = 90. 

6 Conclusion 

A microphone subjected to a plane wave has been studied. Using a propaga- 
tion model in a fluid domain located close to the microphone membrane and 
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measuring the acoustic velocities on the boundary on the domain provides an 
estimate of the acoustic pressure on the microphone membrane. This model 
assumes that the normal acoustic velocity on the microphone membrane is 
uniformly zero and considers the incident wave as plane. 

This preliminary study of the acoustic pressure field in the very close vicinity 
of a microphone highlights an acoustic pressure gradient along the microphone 
membrane in order of magnitude of 0.1 dB. This acoustic pressure difference 
existing between the center and the edge of the membrane can be very prob- 
lematic when microphone calibration requires an uncertainty less than 0.1 dB. 
The hypothesis used in this work (incident plane wave, parabolic profile of 
radial velocity and motionless membrane) tend to minimize this pressure gra- 
dient value. A numerical calculation of integrals (Green fomulation) would 
enable to ignore these hypothesis and could lead to a more realistic pressure 
estimation. 

This preliminary study opens new horizons in microphone calibration research, 
especially in free field conditions. These first results on the acoustic pressure 
field in the very close vicinity set an important question : what is measured by 

the microphone since the acoustic pressure is not uniform along the membrane 

? 

In the future, these results should be validated by measuring the acoustic pres- 
sure field in the very close vicinity of the microphone using, for instance, a 
small probe to not perturb the acoustic field. The proposed model can be im- 
proved by taking into account the membrane motion which could be estimated 
using laser probe measurements. 
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